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0\ '. Abstract 



< 

fSJ ' plasma at zero temperature and high chemical potential to weak color elec- 



Classical transport theory is used to study the response of a non-Abelian 



^ tromagnetic fields. In this article the parallelism between the transport phe- 

ly-^ ' nomena occurring in a non-Abelian plasma at high temperature and high 

' density is stressed. Particularly, it is shown that at high densities it is also 



possible to relate the transport equations to the zero-curvature condition of 
a Chern-Simons theory in three dimensions, even when quarks are not con- 
sidered ultrarelativistic. The induced color current in the cold plasma can 
be expressed as an average over angles, which represent the directions of the 
velocity vectors of quarks having Fermi energy. From this color current it 

X" 

^-H ' is possible to compute n-point gluonic amplitudes, with arbitrary n. It is 

argued that these amplitudes are the same as the ones computed in the high 
chemical potential limit of QCD, that are then called hard dense loops. The 
agreement between the two different formalisms is checked by computing the 
polarization tensor of QED due to finite density effects in the high density 
limit. 
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I. INTRODUCTION 



Quantum Chromodynamics (QCD) undergoes a phase transition at high temperature 
and/or high density [|^. Above a critical temperature and critical chemical potential, quarks 
and gluons are no longer confined. Those extreme regimes are expected to be found in nature 
in some cosmological and astrophysical settings {e.g., in the interior of neutron stars) or in 
heavy ion collisions. These are the regimes of QCD that we are going to discuss here. 

As is well known now, naive perturbative analysis of high temperature QCD fails com- 
pletely. This was realized when physical quantities, such as the gluonic damping rate, were 
found to be gauge dependent when computed following the standard rules of quantum field 
theories at finite temperature. The connection between expanding in loops and expanding 
in the gauge coupling constant g is not valid in this regime. As was realized by Braaten and 
Pisarski 0, as well as by Frenkel and Taylor |Q, there are one- loop corrections, the hard 
thermal loops (HTLs), which are as important as tree amplitudes, and they have to be re- 
summed and included consistently in all computations to non-trivial order in g. Those hard 
thermal loops arise only for soft external momenta (~ gT) and hard internal loop momen- 
tum (~ T), where g <^1, and T denotes the plasma temperature. There is an infinite set of 
HTLs which possess very interesting properties, such as obeying QED-like Ward identities 
and being gauge independent 0, [§]. 

The hard thermal loop resummation techniques that were proposed in were successful 
in providing one-loop gauge independent physical answers. Much related work has been done 
since their discovery (see ref. [||, for a review and references), but here we shall only review 
different approaches to HTLs. 

Taylor and Wong [0 were able to construct an effective action for hard thermal loops 
just by solving the gauge invariance condition imposed on that effective action. Efraty and 
Nair 0, have related that condition to the zero curvature equation of a Chern-Simons 
theory in (2+1) dimensional space at zero temperature, providing therefore a non-thermal 
framework to study a thermal effect. This identification has been used by Jackiw and Nair 
|T0| to derive a non-Abelian generalization of the Kubo formula. 

Other derivations of the effective action of hard thermal loops can be found in the litera- 
ture. Blaizot and lancu |TT| could extract the hard thermal loop effective action by studying 
the truncated Schwinger- Dyson hierarchy of equations, after performing a consistent expan- 
sion in the gauge coupling constant and obtaining quantum kinetic equations for the QCD 
induced color current. Jackiw, Liu and Lucchesi have also shown that HTLs can be 



derived from the Cornwall- Jackiw- Tomboulis composite effective action |]13|, after requiring 
its stationarity. 

A different derivation of the effective action of HTLs has been given in , which doesn't 
make use of quantum field theory, as opposed to all previous approaches. Hard thermal 
effects are exclusively due to thermal fluctuations, and that is why a classical formalism 



to describe them was developed in Just by writing the classical transport equations 
for non-Abelian particles [0, and using an approximation scheme that respects the non- 
Abelian gauge symmetry of the transport equations, the effective action of the infinite set 
of HTLs of QCD could be found. 

A similar situation may be expected to arise for QCD at high density and zero temper- 
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ature [0. Actually, HTLs were also studied when a chemical potential was included for 
quarks, and it was concluded that the only effect of the chemical potential was modifying 
the Debye mass by a term proportional to the chemical potential |T^, |T^. Therefore, at 
very high density or chemical potential fi and zero temperature, one also may expect that 
naive one-loop computations are incomplete, as one-loop diagrams with soft (~ g/j,) external 
momenta, and quarks running inside the loop with Fermi energy, are comparable to tree 
amplitudes, and therefore they would have to be resummed. We shall call those diagrams 
hard dense loops (HDLs) in analogy to the thermal case. 

The purpose of this article is to give a derivation of the hard dense loops of gluonic 
amplitudes of QCD by using the classical transport formalism. We expect that quantum 
field computations of QCD at high chemical potential will reproduce the transport results, 
exactly as happens in the high temperature case, although we do not attempt here to check 
complete agreement between these two different formalisms. 

Although we study only gluonic amplitudes in this article, it should be expected that 
hard dense loops also arise for Feynman diagrams with external quark legs. We do no take 
into consideration those diagrams in the present article, and refer to a recent publication 
T9| which considers the one-loop self-energy of the electron for QED at finite density for 



those kinds of HDLs. 

We have structured this article as follows. In Section II we review briefly the classical 
transport theory for a non-Abelian plasma. In Subsection III A we study how a cold quark 
plasma initially in equilibrium reacts to weak external color electromagnetic fields, and write 
the transport equation that the induced color current density obeys. The parallelism between 
the high temperature and high density cases is clearly established, and we also comment 
about the connection between the transport equation and the Chern-Simons eikonal. In 
Subsection III B we solve the transport equation in momentum space, using recent results 
of ref. In Subsection III C we give the polarization tensor in the quark plasma, and 

prescribe how to derive higher-order point amplitudes. As a check of the agreement between 
the transport approach and the high chemical potential limit of QCD, we compute the 
polarization tensor of QED due to finite density effects in Section IV, and from it we extract 
the HDL corresponding to that graph, which coincides with the one of QCD, up to some 
color factors. Complete agreement is found between the two different computations. On 
this result, we base our expectation that complete agreement for higher n-point gluonic 
amplitudes should also be true. Finally, we present our conclusions. 



II. CLASSICAL TRANSPORT THEORY FOR A NON-ABELIAN PLASMA 



The classical transport theory for the QCD plasma was developed in and further 
studied in |]T^, and here we will briefly review it. Consider a particle bearing a non-Abelian 
SU{N) color charge Q", a = 1, A^^ — 1, traversing a worldline x"{t). The Wong equations 

and (5° (we neglect here the 



(2.1a) 



!1| describe the dynamical evolution of the variables x^. 
effect of spin): 



m ■ 



dr 
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m^^=^7Q"(r)Fr(x)p.(r) , (2.1b) 



m ■ 



dr 
dr 



-gr'^p^{T)A'(x)Q%T) , (2.1c) 



where fate are the structure constants of the group, F^'^ denotes the field strength, g is the 
couphng constant, and we set c = h = ks = 1 henceforth. 

The main difference between the equations of electromagnetism and the Wong equations, 
apart from their intrinsic non-Abehan structure, comes from the fact that color charges pre- 
cess in color space, and therefore they are dynamical variables. Equation ( 2.1c ) guarantees 
that the color current associated to each colored particle, 



j;ix)=gj dTp,{T)Q\T)5^'\x-x{T)) , (2.2) 
where (x^,p^,(5") are solutions of the equations of motion, is covariantly conserved 



{D,3^)^ {x) = d,j^ix) + gUAl{x)rix) = , (2.3) 

therefore preserving the consistency of the theory. 

The usual {x,p) phase-space is thus enlarged to {x,p,Q) by including color degrees of 
freedom for colored particles. Physical constraints are enforced by inserting delta-functions 
in the phase-space volume element dxdPdQ. The momentum measure 

dP = ^^2 9{po)S{p'-m') (2.4) 

guarantees positivity of the energy and on-shell evolution. The color charge measure enforces 
the conservation of the group invariants, e.g., for SU{3), 

dQ = d^Q 6{Q,Q^ - q^) 6{d,tcQ''Q'Q'' - qs) , (2.5) 

where the constants q2 and fix the values of the Casimirs and dabc are the totally symmetric 
group constants. The color charges which now span the phase-space are dependent variables. 
These can be formally related to a set of independent phase-space Darboux variables |]1^ . 
For the sake of simplicity, we will use the standard color charges. 

The one-particle distribution function f{x,p,Q) denotes the probability for finding the 
particle in the state {x,p,Q). In the collisionless case, it evolves in time via a transport 
equation ^ = 0. Using the equations of motion ( |2.1| ), it becomes the Boltzmann equation 

H 



^ -gQ^F^^^^-gU^AlQ^-^ 



fix,p,Q) = 0. (2.6) 



A complete, self-consistent set of non-Abelian Vlasov equations for the distribution func- 
tion and the mean color field is obtained by augmenting the Boltzmann equation with the 
Yang-Mills equations: 
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= J^'^(x) = E ' (2-7) 

species helicities 

where the color current j^"'{x) for each particle species is computed from the corresponding 
distribution function as 

3^''{x)=g J dPdQp>^Q''fix,p,Q) . (2.8) 

Notice that if the particle's trajectory in phase-space would be known exactly, then Eq. 
( p.8|) could be expressed as in Eq. ( p.2|) . Furthermore, the color current (|2.8| ) is covariantly 



conserved, as can be shown by using the Boltzmann equation [0 . 

The Wong equations (|2.1|) are invariant under the finite gauge transformations (in matrix 
notation) 

x'^ = x^ p^=p^, Q = UQU-\ A^ = UA^U-'--U^U-\ (2.9) 

9 oXf, 

where U = U{x) is a group element. 

It can be shown that the Boltzmann equation ( p. 61) is invariant under the above 



gauge transformation if the distribution function behaves as a scalar 

nx,p,Q) = f{x,p,Q) . (2.10) 

To check this statement it is important to note that under a gauge transformation the 
derivatives appearing in the Boltzmann equation (|2.6| ) transform as: 



dx^^ dx^^ \ dx^^ dQ J dp>^ dp>^ dQ dQ 

(2.11) 

that is, they are not gauge invariant by themselves. Only the specific combination of the 
spacial and color derivatives that appears in (|2.6|) is gauge invariant. 

The color current (|2.8| ) transforms under (|2.9|) as a gauge covariant vector: 



f{x) = g j dPdQp^U QU-^ f{x,p,Q) = U f{x)U-^ . (2.12) 

This is due to the gauge invariance of the phase-space measure and to the transformation 
properties of /. 



III. INDUCED COLOR CURRENT IN A COLD QUARK PLASMA 

A. Transport equation for the color current 

In this section we study soft disturbances of a completely degenerate quark plasma. 
We consider a quark plasma at zero temperature and finite density which is initially in 
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equilibrium. In the absence of a net color field, and assuming isotropy and color neutrality, 
the equilibrium distribution function for this system is, up to a normalization constant. 



/W(p) = ^(/i-po) , (3.1) 

where 9 is the step function, and /x is the chemical potential. That is, all particle states are 
occupied with occupancy number one up to the Fermi energy po = A^- 

Now we study how the quark plasma reacts to weak external color electromagnetic fields. 
We consider that the chemical potential is large fi >> fic, where fic is the critical chemical 
potential, while the color fields are soft or of the scale order F^^, ~ gfi'^, where g is the 
coupling constant which is assumed to be small. The distribution function can be expanded 
in powers of g as: 

f = f^'^+9f^'''+9'f^'^ + ..., (3.2) 
The Boltzmann equation ( p.6|) for f^^^ reduces to 



p' {■i-.-sr'-KQc^) f''\-,p^Q) = p'QaF;.^/'\po) . (3.3) 

Notice that a complete linearization of the equation in ^4° would break the gauge invariance 
of the transport equation, which is preserved in this approximation. But notice as well that 
this approximation tells us that f^^^ also carries a (^-dependence. 

We can get the equation that the color current density Jl^{x,p) obeys by multiplying ( p.3|) 
by and Qa and then integrating over the color charges. For quarks in the fundamental 
representation 



/ 



dQQaQb = l5ab. (3.4) 



Taking also into account 

-j-difi-po) = -S{fi-po) : (3.5) 
dpo 

we finally get, after summing over helicities, 

[p-D J^{x,p)r = -g'p''p''F^,{x)6{fi-po) . (3.6) 

The total induced current density can be obtained after summing over the different contri- 
butions due to different quark flavors. In order to simplify the notation we will not write 
the quark flavor index in the color current, chemical potentials, masses, etc., and it should 
be understood that a sum over quark flavors is to be taken in all final formulas. 

From Eq. ( |3.6| ) we see that only quarks which are on the Fermi surface contribute to the 
induced color current in the plasma. Quarks which are inside the Fermi sea are blocked to 
react to the presence of external fields due to the Pauli exclusion principle. Furthermore, 
the plasma only responds to the presence of external color electric fields, and that is why 
only those get a screening mass in this approach. 
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In order to solve ( |3.6| ), we first divide the equation by po. Then we integrate it over pq 



and |p|, using the momentum measure dP ( |2.4| ). Due to the delta function in the momentum 



measure and the delta function in the r.h.s of Eq. ( |3.6| ), J^(x,p) is only non- vanishing when 
Po = v IpP + = A*- Then we can write Eq. ( |3.6| ) as 



[vD J^{x,v)]'' = -M'v^v^F^oix) , (3.7) 

with 

= /fl , (3.8) 
where pp = Vf^'^ ~ is the Fermi momentum and = (1, ^v), and we have defined 

Ja'M = | IPNIp'^^° 2%o)^(p^-m^) J!^{x,p) , (3.9) 
so that the total color current is obtained by integrating over all directions of the unit vector 



J!:ix)= i^j^x^v). (3.10) 



Equation ( p.7|) has the same structure as the transport equation that the induced color 



density quantity J'^^'^\x, V) in a softly perturbed hot quark-gluon plasma obeys |]T^| 

[V ■ D J>'^''\x,V)]'' = -m\V^VP Fl^ix) , (3.11) 

where = g'^T'^ {N + N f / 2) / ?> is the Debye mass squared for a SU{N) non-Abelian group, 
with Nf fiavors of quarks, and is a light-like four vector. 

Equations (|3^ ) and (|3.11|) exhibit a similar structure, although they correspond to two 
different physical situations. In the case of a hot quark-gluon plasma the contribution to 
the induced color current comes from both quarks and gluons which are thermalized at a 
temperature T. This is refiected in the coefficient of the r.h.s. of (|3.11|) , the Debye mass 
squared. For very high temperatures, we expect that quarks in the plasma move very fast, 
and it is a good approximation to consider that their velocities are ultrarelativistic, as long 
as T 3> m, where m is the quark mass, which is then neglected. Gluons travel at light 
velocities, however, without any approximation. These are the reasons why the vector 
in Eq. ( p.ll| ) is a light-like four vector. In the zero temperature quark plasma that we 



are considering here there are no real gluons, and that is why they do not contribute to 
the induced color current. Furthermore, only quarks with Fermi energy contribute to that 
current. For very large chemical potential, it can be a good approximation to neglect quark 
masses if /i ^ m. If quark masses are neglected, then pp = /x, and Eq. (|3.7|) coincides 
with Eq. ( p.ll[ ), except for the the factors m|, and M^. We have decided, however, to keep 



corrections due to quark masses, which is justified if yU > m ^ g^i. Therefore v'^ is not a 
light-like vector in this more general situation where we are considering massive quarks. 
Before solving Eq. (|3.7|) we would like to comment about the connection between this 



transport problem and the Chern-Simons eikonal. Efraty and Nair M, M have shown 
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that the gauge invariance condition for the generating functional of hard thermal loops 
^HTL of QCD can be formally related to the zero curvature condition of a Chern-Simons 
theory in a (2+1) dimensional space-time at zero temperature. Furthermore, it has been 
shown that this condition can be obtained from the transport equation ( |3.11| ) after assuming 
Ja^^K^) = '^^^ same identifications can be done for the cold quark plasma, after 

assuming that the induced color current can be obtained from a generating functional. The 
same steps that were necessary to show that identification for the thermal problem can be 
repeated here with the only main difference that v'^ is not a light-like vector, when quark 
masses are not neglected. 

For completeness we present that identification below. We first define a new current 
density 

J''%x,v) = J''''{x,v) + M^^'A^oix) . (3.12) 
Then JT^" obeys the equation 

[v ■ D J^{x, v)Y = (y ■ A\x)) . (3.13) 

Now we assume that Jj^ can be derived from a generating functional as 

Equation (|3.13| ) then implies that W{A^v) depends only on = —j=vA^^ with / = Pf/ f^-, 
z.e. W{A,v) = W{A+), and = ^-^^^v^. If we now define 

if I \ 1 



2 V " / ; ' " 

— \ Xn -V ■ X , 0_ = , — - 



X- = d ^ ( Xq — ■ yi\ , d_ = —^v-d, (3.15) 



with V = (1, — ^v), and 

F" = ^-iMX, (3.16) 

then, after a Wick rotation to Euclidean space x+ ^ z, x_ — > z, A^ A^ and the 
identification = Eq. ( |3.13| ) is translated into (in matrix notation, and using 

antihermitian generators in the fundamental representation of SU{N)) 

d,a, - d,A, + g[A,, a,] = , (3.17) 

which corresponds to the zero curvature condition of a Chern-Simons theory in (2+1) Eu- 
chdean dimensional space-time, in the gauge Aq = 0. 
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Solutions of Eq. ( |3.17|) are provided by the Chern-Simons eikonal and were studied long 
ago |2^. The generating functional of the induced color current in a cold quark plasma can 
then be expressed in terms of the Chern-Simons eikonal, exactly as in the case of the hot 
quark-gluon plasma. 

In spite of the very interesting connection between these, in principle, unrelated problems, 
we will not pursue the Chern-Simons approach to the transport equations in the remaining 
part of this article. 



B. Solution to the transport equation 



Equation (|3.7|) can be solved after inverting the v ■ D operator and imposing proper 
boundary conditions. Then the solution of (|3.?| ) is expressed in terms of link operators [0 , 
||1 1|| . In this Subsection, we prefer to solve the transport equation by going to momentum 
space exactly as has been done for Eq. (|3.11| ) in ref. pO|- Writing 



Eq. (|3.7|) becomes in momentum space 



d^x 
(2^ 



(3.18) 



v-kj^{k,v)+igfabcj 



V ■ kAl{k) -kov- A'^ik) + igUcj ^ V ■ A\k - g)A^(g) 



(3.19) 



Now, after assuming that J'^{k,v) can be expressed as an infinite power series in the gauge 
field y4^(fc), Eq. (|3.19|) can be solved iteratively for each order in the power series. We will 
impose retarded boundary conditions by the prescription po ^ Po + i^, with e — > 0"*", that 
should be understood in all following formulas. 
The first order solution is 



J,^^'\k,v) = M'v^ (ko ^-^^ - Alik) 



(3.20) 



Inserting ( p.20D in ( p.l9| ) allows solving for the second order term in the series, which reads 



7'^^'Ukv)--taM'f, I v-A\k-q)v.A'^{q) 
The n-th order term (n > 2) can be expressed as a function of the (n — l)-th one as 



(3.21) 



abc 



d^q vA\k-q) 
(27r)4 Vli 



J. 



/t(n-l) 



(3.22) 



The complete expression of the induced color current in the cold quark plasma is thus 
given by 
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(2vr) 



n=l 

'4U oo 



An J (277 ) 



e 



-^k■.Y^J,in)^^^^■^ • (3.23) 



n=l 



The color current is then expressed as an average over the three dimensional unit vector v, 
which represents the directions of the velocity vectors of quarks with Fermi energy. 

It should be possible to construct a generating functional T^dl that generates this 
current by solving 

Jliix) = . (3.24) 

It is first required to check that J^i^) obeys integrability conditions, as has been done in 
ref. [0 for the massless and thermal case. 



C. Polarization Tensor 



The polarization tensor 11^^ can be computed from (|3.20|) by using the relation 

J!:^'Hk) = Kak)Al{k) . (3.25) 

It reads: 



(3.26) 



where we recall that f'^ = (1, ^v). To avoid the poles in the above integrand, we impose 
retarded boundary conditions, i.e., we replace fco by ko + ie. 
The polarization tensor (|3.26|) obeys the Ward identity 



k n'"' 



0, 



(3.27) 



and it is also gauge-independent due to the gauge invariance of the transport formalism. 
The real part of the polarization tensor is 



Reno°(fco,k) = 5,fenKfco,k) 



Re n°Ufco, k) = 5ab ko Mko, k) 



Ren^^b(^o,k) = 6ab 



- ^ n,(A:o, k) + ^ ^ n,(A:o, k) 



(3.28a) 
(3.28b) 
(3.28c) 



where 



9 



\PF2\k\ 



In 



nt(/co,k) 



kP 



I [pf |k| /i /C( 



2 \ fi ko Pf \k\l 



^ko +pF\k\ 



fxko — pF\k\ 



(3.29a) 
(3.29b) 



The imaginary part of the polarization tensor describes damping in the quark plasma, 
explicitly: 



Imn^,^(fco,k) = -6,, M'nkof '^v^v^5{k, - ^k 

J 47r /i 



(3.30) 



Our result of the polarization tensor agrees with the one computed in ref. |2^ using 
classical transport theory for a quark-gluon plasma at finite temperature and baryon density 
in the zero temperature limit. 

Notice that if quark masses are neglected, then pf = fJ., and then ( p.26|) reduces to the 
same polarization tensor as the one found in a hot quark-gluon plasma 
a different screening mass. 

Higher order n-point functions can be found as 



231, 14 , with 



i=n—l 



SJ^ik) 



i=l 



SAll{q^)---SAtlz\iqn-i) 



(3.31) 



A-=0 



where the delta function on the r.h.s. of (|3.31|) accounts for conservation of momentum, and 
permutations of indices and momentum have to be taken. 

Due to the fact that the induced color current can be expressed as an infinite power 
series in the gauge field A'^, it is obvious that there are n-point gluonic amplitudes, with 
arbitrary n, from 2 to oo. All these gluonic amplitudes describe the polarizability properties 
of the non-Abelian plasma. 



IV. THE POLARIZATION TENSOR OF QED AT FINITE DENSITY 

In this Section we present the computation of the polarization tensor of QED due to the 
presence of a background density of electrons, and from it we extract the HDL corresponding 
to that amplitude in the massless case. This computation has been previously performed by 
several authors in Euclidean space-time, see ref. ||2^, p5[, [^], and also [^] and ||29| for 



the non-relativistic case. Here we review how this computation is performed in Minkowski 
space-time when retarded boundary conditions are imposed. We find it more appropriate to 
carry out this computation in Minkowski space-time to obtain the correct analytic properties 
of the retarded polarization tensor. 

In the presence of a finite density of electrons the usual definition of creation and annihi- 
lation operators for particles and "holes" , or antiparticles, of the vacuum has to be modified 
The ground state I'I'o > of the system is constituted by the Fermi sea, that is, by 



electrons occupying all particle states according to the Pauli exclusion principle up to the 
Fermi energy po = fi. Then one defines a creation operator &J^(p) that acting on the ground 
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state creates one particle with energy pq > /i, and the creation operator dl^{—p), for po < /i, 
which creates a "hole" , which may be interpreted as the absence of one electron in the Fermi 
sea. The corresponding particle and hole destruction operators annihilate the ground state 

b^{p) \% >= d^{p) \% >= , (4.1) 

since the ground state doesn't contain electrons with energies above the Fermi energy or 
"holes" inside the Fermi sea. 

The usual averages of products of creation and annihilation operators that are required 
to compute Feynman amplitudes are then modified from their vacuum values, explicitly, 

< $o| ba{p)blip) >= diPO - /U) Sa,a' ■ (4.2) 

Then one can compute the Feynman propagator to get 



zSF{x,y) = < <l>o|r^(x)^(y)|$o > (4.3) 

= e{x^ - y°) / 7^\{i-P + m) e{po - /i) e-'p("-^) 



(27r)3 2po 
J (27r)3 2po L 



where po = y |pp + m^. Notice that in the limit /x — >■ 0, (|4.3| ) agrees with the Feynman 
propagator in the vacuum. 

In order to compute the polarization tensor of QED in the presence of a background 
density of electrons, the propagator ([4.3|) is needed. Retarded boundary conditions are 
taken into account by introducing convergence factors e^^°, e~^^°, with e — >■ 0^, as needed 
when going from configuration to momentum space (see ref. |^ and [1^ for explicit details). 
Taking that into account, we find for the chemical potential dependent part of the retarded 
polarization tensor in the massless case 

d^q 1 1 



(27r)3 2po 2g, 







(^(/x - go) - ^(m - Po)) (4.4) 



- ^(/U - Po) : ; - ^(/i - Qo) 



, . ^yr^ , ■ (4.5) 

Po + go - fco - «e Po + qo + ko + i€_ 

where p = (po,p), p' = (Po,-p), 9 = (go,q), q' = (go,-q), with po = |p|, go = |q|, and 
p = q + k. We have further defined 

T'^'ip, q) = Tr [7'^7 ■ p-f^-f ■ q] . (4.6) 

In the high chemical potential limit, one considers k^, |k| ^ fi. Then one can make the 
approximations @, 

T'"'{p,q) r^Sq^V^'V (4.7) 

T^^'ip, q') ~ T^^'ip', q) ~ Aq^ {v^'V + - 2g>"') (4.8) 

Po- qo- ko -k-V (4.9) 

Po + go ± A;o ~ 2go (4.10) 

e{fi - go) -0{fi-po)^V-k 6{fi - go) , (4.11) 
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where V = (1, g), V = (1, -g) to finally get 

n-(A:) = -9'j^, {V^V^ + V^V^ - 2 ko ^^^^ ^(/^ - ^o) , (4.12) 

which after performing the integral over the modulus of |q| reduces to the same polarization 
tensor ( p.26| ) for the Abelian case and massless case. 



V. CONCLUSIONS 

In this paper we have used classical transport theory to study the response of a completely 
degenerate non- Abelian plasma at zero temperature to weak color electromagnetic fields. We 
have used an approximation scheme for the classical transport equations that respects their 
non- Abelian gauge symmetry, and from it we have obtained the induced color current in the 
non- Abelian plasma. 

The study of the transport phenomena occurring at high density in a quark plasma that 



we have carried out here is quite similar to the one that was done in |TJ] for a QCD plasma 
at high temperature, and we direct the reader to consult that reference. Throughout this 
paper we have stressed the parallelism that these two different problems exhibit, as well as 
their differences. 

Both at high temperature and high density, the color constituents of the non-Abelian 
plasma are not confined, and that is why one can model color degrees of freedom classically. 

The transport equations obeyed by the induced color current for a non-Abelian plasma 
at high temperature and at high density present a similar structure. If quark masses are 
neglected, then the two transport equations are essentially the same, except for the fol- 
lowing. At high temperature T, both quarks and gluons which are thermalized at the 
temperature T contribute to the induced color current in the plasma. At high chemical 
potential, only quarks on the Fermi surface contribute to that current, as there are no real 
gluons at zero temperature. Then, in that situation, the only difference between the two 
transport equations and their solutions comes from the screening masses squared which are 
m% = g^T^{N + Nf /2)/3 and = respectively, for a SU{N) theory with Nf 

fiavors of quarks (and assuming that for all quark fiavors fif = fi). 

We have decided to keep corrections due to quark masses in the present study. The 
solution to the transport equation can be expressed in terms of an average of a three- 
dimensional unit vector v, which represents the directions of the velocity vectors of quarks 
which have Fermi energy. 

We have also shown that the solution of the transport equation can be expressed in terms 
of the Chern-Simons eikonal, even for massive quarks. Even though one cannot define natural 
light-cone coordinates when quarks are not consider ultrarelativistic, we have seen that it 
is still possible to relate the transport equation to the zero-curvature condition of a Chern- 
Simons theory in (2+1) Euclidean space-time. The transport equation only depends on the 
projections of the vector gauge fields over the four vectors v'^ and v'^ defined in Subsection 
III A, and this is enough to define a two dimensional plane where the zero-curvature condition 
of a Chern-Simons theory can be recognized, after proper identifications. 
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We have presented a systematic way of computing n-point gluonic amplitudes using the 
solution of the classical transport equation, and we claim that these results should agree with 
the high density limit of a quantum field theory approach. Complete agreement is found 
for the polarization tensor [n = 2) in the massless case, and complete agreement should be 
expected for higher ?7,-point functions, although we have not check this statement, not are 
we aware of the computation of those amplitudes in the literature. 

Therefore, we claim that a Braaten-Pisarski resummation procedure should also be used 
in QCD with high chemical potential, as was generally expected, although a more detailed 
quantum field theoretical study should be carried out to see exactly how this resummation 
is implemented. We think that it is interesting to have this alternative computation of 
HDLs, specially due to the simplicity and transparent physical interpretation of the classical 
transport formalism. 

Finally, we should mention that although the response theory of the high temperature 
limit (at zero chemical potential) and the high chemical potential limit (at zero temperature) 
of QCD seem almost identical, and one could naively conclude that in both situations static 
color electric fields are exponentially damped (at least at leading order), this is not so. 
Kapusta and Toimela have pointed out that the static potential between two static 



charges in a plasma at zero temperature becomes oscillatory and vanishes as a power at 
large distances due to the existence of a sharp Fermi surface. This effect was known as 



Friedel oscillation in non-relativistic quantum theory |29[ . 

A further study of the high density limit of QCD should be carried out, and we postpone 
that study to a further publication. 
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